function Const=SteadyState_Const(PP,SS)

%% Transition Matrix from V2EV
V2EV        =   sparse(kron(SS.FunApp.ExoState.TrMat',eye(SS.FunApp.V.Num)));

%% Transition Matrix from VCoef to UCoef
f               =   SS.FunApp.U.State(:,1);
Lag_zoID        =   SS.FunApp.U.State(:,2);
N_State         =   length(f);

Lag_oID         =   SS.FunApp.ExoState.State(Lag_zoID,2);
Lag_o           =   PP.ExoState.IdioRet.Node(Lag_oID);
Temp_dE         =   ones(N_State,1)*SS.dE;
Temp_dE(Lag_o==0)=  0;
f               =   f.*(1+Temp_dE)/(1+SS.Pir);
TempBasMat      =   funbas(SS.FunApp.V.Space,f);
[II,JJ,VV]      =   find(TempBasMat);
BasMat_V        =   sparse(II,JJ+(Lag_zoID(II)-1)*SS.FunApp.V.Num,VV,N_State,SS.FunApp.V.Num*SS.FunApp.ExoState.N);
CoefMat_V2U     =   V2EV*BasMat_V;
%% Collect the Results
Const       =   struct('TrMat_V2EV',V2EV,'CoefMat_V2U',CoefMat_V2U);